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We examine superconductivity in layered systems with large Fermi-surface splitting due to coex- 
isting ferromagnetic layers. In particular, the hybrid ruthenate-cuprate compound RuSr2GdCu208 
is examined on the coexistence of the superconductivity and the ferromagnetism, which has been 
observed recently. We calculate critical fields of the superconductivity taking into account the Fulde- 
Ferrell-Larkin-Ovchinnikov state in a model with Fermi-surfaces which shapes are similar to those 
obtained by a band calculation. It is shown that the critical field is enhanced remarkably due to 
a Fermi-surface effect, and can be high enough to make the coexistence possible in a microscopic 
scale. We also clarify the direction of the spatial oscillation of the order parameter, which may be 
observed by scanning tunneling microscope experiments. 



Recently, coexistence of superconductivity and ferro- 
magnetism has been reported in the hybrid ruthenate- 
cuprate compounds i?i.4Ceo.6R.uSr2Cu20io-<5 (R = Eu 
and Gd) and RuSr 2 GdCu2 0s These compounds 

have similar crystal structures to the high-T c cuprate 
superconductor YBa2Cu307 except that layers of CuO 
chains are replaced with ruthcnatc layers. Experimen- 
tal and theoretical studies indicate that the ruthenate 
layers are responsible for the ferromagnetic long range 
order [^4j] , while the cuprate layers for the superconduc- 

One of the remarkable features of these compounds is 
that the superconducting transition occurs at a temper- 
ature well below the ferromagnetic transition tempera- 
ture unlike most of the other ferromagnetic supercon- 
ductors. For example, in RuS^GdC^Og, the supercon- 
ducting transition was observed at T c ~ 46K, whereas 
the ferromagnetic transition at Tm ~ 132K There- 
fore, the ferromagnetic order can be regarded as a rigid 
back ground which is not modified very much by the ap- 
pearance of the superconductivity. This picture is also 
supported by experimental observations Jl] || . 

According to the first principle calculations by Pickett 
et al. [Q , magnetic fields in the cuprate layers due to the 
ordered spin moment in the ruthenate layers are much 
smaller than exchange fields mediated by electrons. The 
exchange fields play a role like magnetic fields which act 
only on the spin digrees of freedom but do not create 
Lorentz force. Therefore, the present system is approx- 
imately equivalent to a quasi-two-dimensional system in 
magnetic fields nearly parallel to the layers. 

However, such Fermi-surface splitting gives rise to pair- 
breaking effect as well as that due to a parallel mag- 
netic field. The exchange field in RuS^GdC^Os is very 
large and seems to exceed the Pauli paramagnetic limit 
(Chandrasekar-Clogston limit) ||. The Pauli paramag- 
netic limit Hp at T = is roughly estimated from the 
zero field transition temperature Tc°^ by a simplified for- 
mula fi e Hp = 1.25 T c , where fi e denotes the electron 
magnetic moment. For RuS^GdCt^Os, since the ex- 



change field exists in practice, T c of isolated cuprate 
layers is not known, but it will be appropriate to as- 
sume Tc°^ < 90K from the transition temperature of 
YBa2Cu307+5 at the optimum electron density. Hence 
we obtain /i e Hp < 110K at T = from the above for- 
mula. On the other hand, the band calculation gives an 
estimation fi e B cx = A cx /2 ~ 25meV/2 ~ 107K 0. It is 
remarkable that the superconducting transition occurs at 
such a high temperature T c s» 46K in spite of the strong 
exchange field of the order of the Pauli paramagnetic 
limit at T = 0. 

There are some mechanisms by which the critical field 
of superconductivity exceeds the Pauli limit. For exam- 
ple, the triplet pairing superconductivity is an impor- 
tant candidate. However, from their crystal structures 
and high transition temperatures, it is plausible that the 
present compounds are categorized as high-T c cuprate su- 
perconductors and therefore the superconductivity is due 
to an anisotropic singlet pairing with line nodes, which 
is conventionally called a d-wave pairing. For the singlet 
pairing, possibility of an inhomogeneous superconducting 
state that is called a Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO or LOFF) state was discussed by Pickett et 
al. as a candidate for the mechanism. 

On the possibility of the FFLO state, they pointed 
out that there are nearly flat areas in the Fermi-surfaces 
in RuSr 2 GdCu 2 8 , which favor the FFLO state. It is 
known that the FFLO critical field diverges at T = in 
one dimensional models. However, if the Fermi-surfaces 
are too flat, nesting instabilities, such as those to spin 
density wave (SDW) and charge density wave (CDW), 
are favored for realistic interaction strengths. For the 
present compound, the nearly flat areas are not so flat 
that the nesting instabilities occur, but the small curva- 
ture still enhances the FFLO state ||. 

It is also known that even in the absence of the flat ar- 
eas, the critical field is enhanced in the two-dimensional 
(2D) systems in comparison to the three dimensional sys- 
tems P,[s| JlC]jri"| . Further, when the Fermi-surface struc- 
ture of the system satisfies a certain condition, the critical 
field can reach several times the Pauli limit even in the 
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absence of nearly flat areas |U2[ . Such a Fermi-surface ef- 
fect can be regarded as a kind of nesting effects analogous 
to those for SDW and CDW §. The "nesting" effect 
was examined in details in our previous papers, where 
2D tight binding models are studied as examples ■ 

Direct evidence of the FFLO state may be obtained by 
scanning tunneling microscope (STM) experiments. For 
a comparison with experimental results, spatial structure 
of the order parameter should be predicted theoretically. 
In particular, direction of the modulation of the order 
parameter is important. It may appear that the modula- 
tion must be in the direction perpendicular to the flattest 
area of the Fermi-surface, because then the spatial vari- 
ation is minimized. However, in some of 2D models, it is 
not perpendicular to flattest areas jn]|^] . Only explicite 
calculations which take into account the Fermi-surface 
structure could clarify the direction of the modulation. 

Therefore, the purposes of this paper are (1) estima- 
tion of the critical field of superconductivity including 
the FFLO state to examine the possibility of coexistence 
of singlet superconductivity and ferromagnetism in a mi- 
croscopic scale, and (2) clarification of the direction of 
the spatial oscillation of the order parameter to compare 
with results of STM experiments possible in the future. 
We examine a tight binding model with Fermi-surfaces 
which shapes are similar to those of RuSr2GdCu2 0s, be- 
cause the quantities that we are calculating are sensitive 
to the Fermi-surface structure. 

Recently, the FFLO state has been studied in a 
tight binding model with only nearest neighbor hop- 
ping (1^|l| ■ It was found that ratio of the FFLO critical 
field and the Pauli limit is small near the half filling. Zhu 
et al. have discussed that hence the coexistence of the 
superconductivity and the ferromagnetic order is diffi- 
cult except in the vicinity of the ferromagnetic domains 
near the half filling ||@. However, some experimen- 
tal results indicate coexistence in a microscopic scale and 
a bulk Meissner-state Here, we should note that 

the tight binding model with only nearest neighbor hop- 
ping can not reproduce the shapes of the Fermi-surfaces 
of RuSi^GdC^Os. By taking into account the realistic 
Fermi-surface structure, we will show below that the crit- 
ical field is enhanced remarkably and thus the coexistence 
in a microscopic scale is possible in this compound. 

First, we define the tight binding model 



Ho = V e 
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with a dispersion relation 

e pcr = — 2t(cos p x + cos p y ) — 4^2 cosp x cosp y — n+ha, (2) 

where h denotes the exchange field. When we apply the 
present theory to type II superconductors in a magnetic 
field B, h is written as h = fi e \B\. We use a unit with 
t = 1 and the lattice constant a — 1 in this paper. 

We take the value of the second nearest neighbor hop- 
ping energy t 2 — — 0.6i, which gives shapes of the Fermi- 



surfaces similar to the symmetric Cu02 barrel Fermi- 
surfaces obtained by Pickett et al. |Q] at n = 1.1 as shown 
in Fig.[l| Here, n is the electron number par a site. 
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FIG. 1. Fermi-surface of the present model Hamiltonian for 
t2 = —0.6 and n — 1.1. 

We calculate the critical field in the ground state for 
n = 0.92 ~ 2, applying a formula developed in our pre- 
|l2|l ■ For anisotropic pairing 



vious papers [[PL 

A(p,r) = A Q7Q (p)e^ r 
(p = p/|p|), the critical field is give by 



(3) 



h c = max 
q 



■ exp 



d P\\ Pl(0,P\ 
2tt N a (0) 



log 1 



2/i c 

(4) 

where A Q o = 2wd exp(— l/g a N a (0)) s» 1.76&bT c and 
p±(0,p\\) = /9_l(0,P||)[7 q (p)] 2 with the momentum de- 
pendent density of states pj_(e,pn). Here, pii denotes 
the momentum component along the Fermi-surface. The 
pairing interaction is assumed to have a form 

V(p,pO = -ffa7a(p)7o(p')' ( 5 ) 

In particular, for d-wave pairing, we use a model with 

7d(p) oc cosp x - cospy, (6) 

where p x and p y are the momentum components on 
the Fermi-surface in the directions of p. In our previ- 
ous papers, it was shown that the qualitative and semi- 
quantitative results are not sensitive to details of the form 
of 7d(p) 1 1 1 1^] . An effective density of states N a (0) for 
anisotropic pairing is defined by 



N a (0) = N(0)([ la (p)} 2 
with an average on the Fermi-surface 
dp\\ P±(Q,P\\) 
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where N(0) is the density of states at the Fermi level. 
The Pauli limit Hp for anisotropic pairing is calculated 
by 



with 
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exp ( 



y/([ 7tt (p)] 2 ) A a0 
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[ 7a (p)] 2 log[l/| 7a (p)|]) , 
(Mp)] 2 ) 
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In the above equations, the vector q is the center-of- 
mass momentum of Cooper pairs of the FFLO state. 
From the symmetry of the system, there are four or 
eight equivalent optimum vectors (q m 's), depending on 
whether q is in a symmetry direction or not, respectively. 
Actually, arbitrary linear combination of exp(iq m • r) 
gives the same second order critical field, and the de- 
generacy is removed by the nonlinear term of the gap 
equation below the critical field [|7|Jl6|]. However, regard- 
ing the critical field and the optimum direction of the 
oscillation of the order parameter near the critical field, 
it is sufficient to take a single q as in eq.(3). 

Figures || and |^ show numerical results of the critical 
fields for t<i = —0.6, with our previous results for t% = 
(dotted lines) |f^]. It is found that the critical fields are 
remarkably enhanced near the electron densities n « 1.46 
and 1.20 for the s-wave and the d-wave pairing, respec- 
tively. For example, at the electron density n — 1.1, the 
ratios of the critical field to the Pauli paramagnetic limit 
are approximately equal to 1.66 and 3.19 for the s-wave 
and the d-wave pairing, respectively. These values (es- 
pecially the latter) seem to be large enough to make the 
coexistence possible in RuS^GdC^Os- 

In Fig for the d-wave pairing, both the critical fields 
for <p q = 7r/4 and <pq = are shown, but the highest one 
is the final result of the critical field given by eq. (4) . Here, 
</? q is the angle between the optimum q and one of the 
crystal axes. It is shown by a numerical calculation that 
the critical fields for the other values of ipq are lower than 
the higher one of the critical fields for ipq — ir/4 and 0. 
Thus, the direction of the optimum wave vector q jumps 
from ifq = 7r/4 to ifq — at n w 1.63. On the other hand, 
for the s-wave pairing, tpq = 7r/4 is the optimum in the 
whole region of the electron density. These behaviors are 
different from that for t% = 0, in which tpa = [I4 ]. 

For t2 — —0.6, a cusp is seen in Fig.[2| for the s-wave 
pairing, whereas it does not appear in Fig.|3| for the d- 
wave pairing. The physical origin of the cusp at n w 
1.46 is that the Fermi surfaces satisfy a certain condition 
there, which was explained in our previous paper for £2 = 
. It is related to how the two Fermi-surfaces touch 
by the translation by the optimum q. In the present 
case (t2 = —0.6 and n ~ 1-46), the touch occurs in the 
(110) direction, but because of the nodes of the order 
parameter the "nesting" is not efficient for the d-wave 
pairing. Therefore, cusp does not appear for the d-wave 
pairing. 




FIG. 2. Critical fields of the FFLO state of the s-wave 
pairing for n = 0.92 ~ 2 at T — 0. Solid and broken lines 
show the results for ti = —0.6 and ti = 0, respectively. 
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FIG. 3. Critical fields of the FFLO state of the d-wave 
pairing for n = 0.92 ~ 2 at T — 0. Solid and broken lines 
show the results for <p q = 7r/4 and ipq — 0, respectively, when 
t2 — —0.6. The dotted line shows the result for ti = 0. Thin 
solid line and thin dotted line show the Pauli paramagnetic 
limits in the unit of A^o, for = —0.6 and 0, respectively. 

In spite of the absence of cusp behavior, the critical 
field is still very large for the d-wave pairing near the 
half-filling. Figure |] shows the nesting behavior of the 
Fermi-surfaces at ti = —0.6 and n = 1.1. The direction 
of the optimum vector q is (^ q = tt/4, and the Fermi- 
surfaces touch at two points (i.e., two lines in the p x PyPz~ 
space), (p x ,p y ) w (1.113-7T, 1.71371") and (1.7137T, 1.1137r). 
Since y q = 7r/4 is also the direction of a node of the d- 
wave order parameter, it may appear that this direction 
is less favorable. However, in actuality the critical field 
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is remarkably enhanced for this "nesting" vector q, since 
it gives two nesting lines which are far away from the 
nodes but near the flattest areas, as shown in Fig.^. Be- 
sides, they are near both the maxima of the d-wave order 
parameter and the van Hove singularities, which also en- 
hance the critical field. As the electron density increases, 
the two nesting lines approach to the line node of the 
order parameter, and thus the critical field decreases. 
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FIG. 4. Fermi-surface nesting for the FFLO state and the 
optimum wave vector q of the FFLO state at the critical field 
(h c = 1.87Ado) at n = 1.1. Solid and broken lines show the 
Fermi-surfaces of the up and down spin electrons, respectively. 
Dotted line shows the Fermi-surface of the down spin electrons 
shifted by q. Small arrows show the wave vector q. They 
are placed at the momenta at which the two Fermi-surfaces 
touch. Dotted broken lines show the nodes of the d-wave 
order parameter. We use a large value of A^o = 0.3t/1.87 
(i.e., h c = 0.3t) in order to make the displacement visible. 

Since the optimum direction </? q = 7r/4 is in a sym- 
metry line, there are four equivalent directions, that is, 
ifq = ±7r/4 and ±37r/4. Therefore, symmetric linear 
combinations such as 



A(p, r) oc cos(gx') 

A(p, r) oc cos(gir') + cos(gy') 



(11) 



are convincing candidates, which may be observed in the 
present compound, where x' = (x + y)j\f2 and y' — 
(x — y)/y/2. In particular, the 2D structures such as the 
latter of eq.(ll) are favored at high fields Jig] . 

For the FFLO state to appear, temperature needs 
to be lower than the tri-critical temperature T* of the 
FFLO, BCS and normal states. T* is generally equal 
to about 0.56Tc°^ in simplified models such as eq.(5). 
If we apply this to the present system RuSr 2 GdCu 2 08, 

T* > T c « 46K requires T c (0) > 82K. This condition for 

Tc (0) 

may be relaxed by taking into account a mixing of 



order parameters of different symmetries, which increases 
T* f7|. 

In conclusion, the FFLO critical field of the cuprate 
layers is remarkably enhanced by an effect of the Fermi- 
surface structure. The direction of the spatial oscillation 
of the order parameter is in the (110) direction both for 
the s-wave pairing and the d-wave pairing. Although we 
examined only the ground state in this paper, the result 
H c /Hp w 3.19 at T — is large enough to support co- 
existence of the superconductivity and the ferromagnetic 
order in a microscopic scale in RuSr 2 GdCu 2 0s. Calcu- 
lation for finite temperatures is now in progress. 
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